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ABSTRACT 



It is now routine to measure the weak gravitational lensing shear signal from the mean el- 
lipticity of distant galaxies. However, conversion between ellipticity and shear assumes local 
linearity of the lensing potential (i.e. that the spatial derivatives of the shear are small), and this 
condition is not satisfied in some of the most interesting regions of the sky. We extend a deriva- 
tion of lensing equations to include higher order terms, and assess the level of biases introduced 
by assuming that first-order weak lensing theory holds in a relatively strong shear regime. We 
find that, even in a worst-case scenario, a fully linear analysis is accurate to within 1% outside 
~ 1.07 times the Einstein radius of a lens, by deriving an analytic function that can be used 
to estimate the applicability of any first-order analysis. The effect is too small to explain the 
discrepancy between weak- and strong-lensing estimates of the mass of the bullet cluster, and 
should not impact cluster surveys for the forseeable future. In fact, it means that arclets can be 
used to measure shears closer to a cluster core than has been generally appreciated. However, 
at the level of accuracy demanded by future lensing surveys, this bias ought to be considered in 
measurements of the inner slope of cluster mass distributions and the small-scale end of the mass 
power spectrum. Both of these are central in determining the relationship between baryonic and 
dark matter. 

Subject headings: gravitational lensing 
1. Introduction 



Gravitational lensing is the deflection of light 
rays from a background li ght source by an inter- 
venin g gravitational field (|Melherlll999t iRefregierl 
2003). It is one of the most promising probes of 
the distribution of dark matter, and hence the ef- 
fects of dark energy. Along lines of sight where 
the deflection is sufficient, "strong lensing" visibly 
distorts (and often multiply images) the shapes 
of individual background galaxies. However, only 
"weak lensing" is produced along most lines of 
sight, even those passing through the outskirts of 
galaxy clusters. This weaker but ubiquitous signal 
has to be collected statistically. To first order in 
a Taylor series, it is obtained from the mean ellip- 
ticity of an otherwise uncorrela ted set of galaxies 
( Bartelmann fc Schneider 2000h . 

Weak lensing measurements have n ow been well 
used to map the distribution of mass ( Clowe et al 
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that dominate the contri bution to the powe r spec- 
trum on small scales (e.g. Smith et al.ll2003h . Near 
such regions, the first-order assumptions implicit 
in a weak lensing analysis no longer necessarily 
hold. In this paper, we expand the Taylor se- 
ries of the weak lensing equation to include the 
next-highest terms, and investigate the level of 
bias in shear measurements that rely upon simple 
measurements of ellipticity. 

We derive the lensing equations in !j2] We check 
our results using raytraced simulations in Sj3j and 
we discuss their implications in j|4l 
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2. Lensing Transformations 

2.1. The Usual First-Order Treatment 

A general gravitational lens deflects a light from 
a position x' in a background (source) image to a 
position x in the observed (lens) plane, such that 



a[x) , 



(1) 



with a deflection angle predicted by General Rel- 
ativity in the weak field limit of 



a(x) = W(x) 



(2) 



and where \& (a?) is the Newtonian potential of the 
lens, $(x, z), projected onto the plane of the sky. 

Crucially, the gravitational field and the deflec- 
tion angle vary across the sky. Assuming (the lo- 
cal linearity condition) that the change is linear 
on scales the size of a galaxy, it can be described 
to first order by a coordinate transformation 



(3) 



The first derivative term represents an unmeasur- 
able centroid shift. Placing the origin of the co- 
ordinate system at the galaxy's observed center of 
light, we are left with 

x' i = A ij x j + ... , (4) 

where the Jacobian of the transformation is 
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We have introduced the usual notation of conver- 
gence k(x) = V 2 ^(x)/2, which is proportional to 
the distribution of mass projected along a line of 
sight, and two components of shear ji(x). The 
inverse mapping is simply 



(A) 



V 



(7) 



It is always possible to adopt an arbitrary 
choice of rotation for the coordinate system such 
that 72 = (so A is diagonal), and invoke par- 
ity symmetry to consider only that the potential 
increases to the right (hence 71 < 0). We also 



work only in the "positive parity" lensing regime 
(outside the critical curve), where det^4 > 0. Our 
analysis is equally valid inside the critical curve, 
although breaks down if a part of the image c rosses 
the critical curve (c./. ISchneider fe Enl2007l) . 

The shape of a galaxy image I(x') can be quan- 
tified via its intrinsic ellipticity 
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where its quadrupole moments are 
/ I(x') x[ x' 3 d 2 x' 
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Under the (locally linear) lensing transforma- 
tion (|7|), the galaxy's observed ellipticity becomes 



X 



obs 



X j nt + 2 7i - X j nt ( X f7 J ) 



(10) 



to first order in 7 (|Seitz fc Schneideijll995l ). Av- 
eraging over a population of galaxies with uncor- 
rected intrinsic shapes (x mt ) = 0, an estimator 7 
can then recover the gravitational shear signal 



,-v_ (X° bs ) _, 
W - 2 _ ((x i„ t)2) Hi 



(11) 



The variance in the denominator can be closely 
approximated by the observed value. It i s typically 
of order 0.4 ( e.g. iLeauthaud et al~ll2007t) . 

For practical purposes, a weight function W{x) 
with finite support is also usually applied to the in- 
tegrals in equation ([9]). This complicates the shear 
estimator: the shear polarizability tensor P 1 in 
Kaiser. Squires fc Broadhurstl ( 19951 ). which gen- 
eralizes the denominator of equation pip , involves 
derivatives of W{x). However, P 7 is typically fit- 
ted from a large ensemble of galaxy shapes to re- 
duce noise, and almost all of those galaxies will be 
on lines of sight unaffected by higher order lensing 
terms. We therefore ignore the effect hered. 



3 As pointed out during the d erivat ion of "reduced shear" 
by iBartelmann & Schneider] l|200Cir i . a galaxy's flux I(x') 
could be replaced in eq. (0 and throughout by a monotonic 
function of intensity f(I(x')), without any change in the 
formalism. This approximates a useful weighting scheme. 
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2.2. Higher order terms 



Continuing th e Taylor series in equation pt , we 
can write (c.f. iGoldberg fc Nataraian 2002 
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Repeatedly substituting the simple form 
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into itself then yields the inverse mapping 

Xi = (A) 

+ x l2{A)^{A)^{A)i q 1 ^, jM x' p x' q 

+ 1 /s(A)il(A)^(A)Y q 1 (A)^, jHm x' p x' q x' r 

+ 7 2 C^ 1 (A)T* {A)-l {A)-, 1 

The various terms are listed in order of decreasing 
importance. Third de rivatives of \fr are related 
to the flexion si gnal ( Goldberg fc Bacon 20051 : 
Bacon et al.l 12006). This is small for realistic po- 
tentials; higher derivatives of ^ will be smaller 
still. Note that this relation still holds locally even 
if there are multiple images, but that there will be 
different values of A at each image. 

To complicate matters, this mapping now shifts 
the galaxy's center of light. If the background im- 
age were correctly centroided (i.e. (x') = 0), the 
observed centroid would be 
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plus smaller contributions coming from the galaxy's 
intrinsic octopole moment. In a coordinate system 
centered on the observed image, the mapping (for 
a fully general potential) is therefore (c.f. eq. [7} 



In practice, a galaxy's intrinsic quadrupole mo- 
ments cannot be observed. We expand them as 
a function of the galaxy's observed shape using 
equation (TT2"]) . However, several non- negligible co- 
efficients produce an unwieldly general expression. 

To make the equations more tractable, we now 
fix various properties of the lens and the source 
galaxy. We first set to zero all derivatives of \I/ 
that are "odd" at 90° (*,ii 2 , *,222, *,m2 and 
^,1222). For a circular (or nearly circular) poten- 
tial that has been rotated so that 'J, 12 = 0, this 
assumption will be (nearly) accurate. It is also 
explicitly true at the major and minor axes of an 
elliptical potential. 

Since we are in a fairly strong lensing regime, 
it is not unreasonable to assume that 7 ;§> x mt j 
so the galaxy can be considered intrinsically cir- 
cular. It still has a size R 2 = 2Q 1 ^ = 2Qf 2 t and 
concentration index 

jm\x\*d*x 

'~ (R 2 ) 2 jl(x) d 2 f { ' 

which would be 2 for a Gaussian, 10/3 for an ex- 
ponential, and higher still for a de Vaucouleurs 
profile. The observed ellipticity becomes 

a 2 d 2 R 2 



x° hs 



4(a 2 



d 2 ) 2 

:115a 4 *, 2 !! -(12a 2 d 2 + 4ad 3 - Sd 4 )^ 2 122 
-2a 2 rf(2a-3d)*, m *,i22 
+4a 3 *,mi -4ad(a - rf)*,ii22 -4d 3 *, 22 22 } 



} (18) 



where a = 04 _1 )ii = (1 - *,n ) _1 and d = 
(A~ x )22 = (1 — ^,22 ) _1 ar e unitless. For a Singu- 
lar Isothermal Sphere (SIS) lens, ^>(x) — 9e\x\, 

cR 2 [l2r 3 - (7 - \)r 2 - 12r - 7] 
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4(r-l) 4 (r 2 + (r-l) 2 ) 



(19) 

where r = \x\/9e- The deviation from an elliptic- 
ity assuming local linearity, x hn > tends as R 2 /9 2 E . 

3. Verification through raytracing 



(A? *5 ( 16 ) 

l h(A)^(A)^(A)i q H, jkl (x' p x' q - Qj£) 
Va (A)J (A)^ (A)^ 1 (A)£% jMm x' p x' q x' r 
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We have developed a simple raytracing rou- 
tine to deflect rays via equation deforming 
the intrinsic shapes of source galaxies into arcs. 
The upper panel of figure Q] demonstrates the ef- 
fect of a singular isothermal sphere lens with Ein- 
stein radius 9e on an intrinsically circular Gaus- 
sian source with a = 0.016>£. Note that this 
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transformation. 



Discussion 




Fig. 1. — Upper images: The observed shape of an 
intrinsically circular galaxy with a Gaussian radial 
profile and size a = O.OWe, at various positions 
behind a singular isothermal sphere lens. The im- 
ages are presented with a logarithmic color stretch. 
Main panel: The solid lines show the object's ellip- 
ticity predicted by the usual linear model and our 
higher order model. The dotted line shows mea- 
sured values from a fully raytraced simulation. 



is a worst-case scenario in several respects, with 
more concentrated or smaller galaxies being less 
affected. If the lens were Abell 1689, this would 
correspond to a z = 1 g alaxy of FWHM ~ 1" 
(jClowe fc Schneider 2001), wh i ch is amongst the 
largest 1 % of Leauthaud et al.l ( 2007t )'s catalog at 
magnitude i' — 25. 

The main panel of figure [T] shows the measured 
ellipticity of the raytraced images, and the predic- 
tion of linear and higher order models. These con- 
verge away from the lens; the slight difference be- 
tween them and the raytraced version is an effect 
of image pixellization. converging slowly. Near 
the lens, our nonlinear model (TIB]) again presents 
a worst case of the deviation from a linear predic- 
tion. It differs from the raytraced measurements 
due to even higher order terms in the coordinate 



We have derived the next-highest terms in the 
coordinate transformation relevant for weak grav- 
itational lensing, by dropping the assumption of 
"local linearity" , which acts as a useful constraint 
on the applicability of the linear approximation. 
The resulting equations are not elegant, but can be 
simplified by making several reasonable assump- 
tions about the galaxy's intrinsic shape and the 
lens profile. As expected, the perturbations from 
linear lensing theory are greatest for large galax- 
ies; they increase as the size of the galaxy squared. 
Like with gravitational flexion, this is simply due 
to the accumulating change in shear signal across 
the width of an image. 

A linear lensing analysis systematically overes- 
timates the shear signal near the core of galaxy 
clusters. However, even in the worst case scenario, 
it is acceptable surprisingly far into the non- linear 
regime. Assuming a value of 1.6 for the denomi- 
nator of equation (fTTj) . it is within 1% of the true 
shear outside ~ 1.07#£, where 7 ~ 0.47, and the 
reduced shear, g ~ 0.93. Compared to other po- 
tential errors, we therefore conclude that this will 
be of only minor concern for measurements of the 
mass of individual (or even stacked) clusters in 
immediately forthcoming surveys. For example, 
the effect is in the right direction but an order of 
magnitude too small to explain the discrepancy 
between measurements of the mass in the bullet 
cluster ( Clowe et al. 2006) via strong and weak 
lensing. However, the effect ought to be considered 
by programs measuring the inner slopes of cluster 
mass distributions or the mass power spectrum on 
small scales. The effect can become relevant at 
about the level of statistical accuracy proposed for 
next-generation surveys. 

We have not investigated the correction for a 
point spread function or the use of a weight func- 
tion while measuring galaxy shapes. A full analy- 
sis of these would be interesting in future work. 

The authors thank Douglas Clowe, Yannick 
Mellier and Barnaby Rowe for useful discussions. 
This work was supported by NASA grant ATP04- 
0000-0067 and DoE grant FG02-04ER41316. 
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